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Abstract 

> 

'nI" . We investigate the production of the (j) meson and the baryon which interact 

weakly with hot hadronic matter, and thus their spectra reflect the early stage of the 
' heavy ion collisions. Our analysis shows that the hadronization temperature, T/j^^, 

■ and the transverse flow, Vj^, of the initial deconfined phase are strongly correlated: 



ON 



Thad + a • (4)^ = 0.25 GeV, where a = 0.37 GeV in the Pb+Pb colhsion at 158 
^ AGeV/c. When choosing appropriate initial values of T^ad and Vj^ from the temper- 

p • ature region T^ad = 175 ± 15 MeV, the measured p meson spectra was reproduced 

^ , surprisingly well by the MICOR model. We have found weak influence of final state 

Q ' hadronic interactions on the transverse hadron spectra at rriT — rrii > O.SGeV. 

^. ■ 

Recent analysis of NA50 Collaboration on and p/u meson transverse momentum 
I spectra in Pb+Pb collision at the CERN SPS, yielded remarkably small values for the 

\-{ I inverse slopes of these particles, namely T{(f)) = 222 ± 6 MeV and T{p/uj) = 219 ± 5 MeV 

0. Furthermore, the analysis by the WA97 Collaboration 0| reveals a slightly larger 
value for the inverse slope of the f2~ + f2 particle spectra, T{Q) = 238 ± 17 MeV ||^. 

These inverse slopes do not fit the suggested straight line of Tj = Tfo + mi{vT)'^ [|, P], 
which may indicate a mutual hadronic transverse flow and which holds for vr^, , p, p, 
(f), A, A and d data of NA44 [|] and NA49 0, |, H with freeze-out temperature Tfo = 175 
MeV and average transverse flow (vt) = 0.35c. On the other hand, the measured slopes 
of the S~ and did not follow this linear dependence on hadron masses rrii This 
phenomena was explained by the early freeze-out of the S and the Q in an analysis based 
on the RQMD model [l^, where final state interactions were considered to reproduce 
the measured hadronic transverse slopes. Since the interaction of the meson with non- 



strange baryon-rich matter is also small [0, the measured value of T((f)) from NA50 
agrees with this picture. 

However, accepting that the momentum spectra of the and Q reflect the early stage 
of the heavy ion collisions, we can ask if this early stage could have been a deconflned 



matter, or simply resonance matter as indicated by RQMD [|T0 . 



In this paper we investigate the direct production of the Q, and p from thermahzed 
massive quark matter. We use the Microscopical COalescence Rehadronization (MICOR) 



model |]T2|, |13l O], which is the successor of the ALCOR |T5|, |T6l and the Transchemistry 



models fl^- The ALCOR and Transchemistry models were constructed to determine 
the total number of hadrons produced from quark matter. The MICOR model is able 
to determine the full momentum spectra of final state hadrons. For the investigation 
of the production of (p and Q we use a part of the full MICOR model, which part will 
be summarized here. A full scale description of the MICOR model and the detailed 
investigation of all particle spectra will be published elsewhere [0. 

One description of the deconfined phase consists of massive quarks, antiquarks and 
gluons in the temperature region 1 < T/T^ < 3 [jl9|. In the MICOR model, which is 
based on this description, the massive gluons (whose number is suppressed relative to 
the quarks and antiquarks) become responsible for an attractive effective potential be- 
tween the colorful quarks and anti-quarks and this effective strong interaction drives the 
hadronization via the coalescence of quarks. The obtained quark-antiquark plasma is as- 
sumed to be a fully thermalized state, and in the hadronization stage it is characterized by 
a hadronization temperature Thad, a transverse flow f j. and a Bjorken-scaled longitudinal 
flow. In spite of the thermal initial conditions, the coalescence of the massive quarks and 
antiquarks produces a hadron resonance gas which is out of equilibrium. Furthermore, 
the decays of the resonances do not lead to equilibrium, but keep the final state hadrons 
in their out of equilibrium state. While the final state hadronic spectra can be fitted by 
inverse slopes familiar from equilibrium descriptions, this does not mean that we have 
equilibrated hadronic final states in our calculation. 

During hadronization meson-like objects are formed in one step via quark-antiquark 
coalescence. Baryon-like objects are formed in two steps: the formation of diquarks 
from quarks are followed by the coalescence of diquarks and quarks into baryons. The 
presence of strongly correlated hadron-like objects in the deconfined phase is supported by 



lattice-QCD results |20]. We use the name 'prehadron' for such a correlator. The values 
of prehadron masses, rripreh, are not sharply determined. In the static case one may 
determine the spectral function of the correlator, extracting some 'mass' and 'width'. 
However, in our case the two-body coalescence process dynamically creates an off-shell 
prehadron with acoai cross section (see Ref. [|T5| for the coalescence cross section). 

We assume that during the phase transition period these color-neutral prehadrons 
can escape the deconfined state without disintegration and they will become the species 
of the produced excited hadronic gas. Large constituent quark masses imply a further 
assumption that the production of the excited J=l vector mesons {p,K*,(j)) and J=3/2 
baryons (A, S*, S*, Q) are favored during the hadronization. The J=0 pseudoscalar 
mesons (pions and Kaons) can not be formed by quark coalescence, but appear together 
with the J=l/2 baryons as products of the decay of the heavier excited hadrons. The 
expected minor difference between the prehadron masses and the excited hadron masses 
will be corrected in a last step: when the prehadron escapes from the deconfined region, 
it becomes on-shell by conserving its velocity. 

In general, the coalescence-type hadron production, qi + q2 — > h, can be described 
by a relativistic rate equation based on the densities pl[]: 



'?1>92 



Here n^^ and the local quark densities, Uh is the prehadron density, is the 

four-velocity of the matter, f 12 is the relative velocity of the two quarks, cTg^g^ the quark 
coalescence cross section and {<Jg^g^Vi2) is the momentum space average of their product: 



For the invariant quark distributions, fi{pi,Xi), thermalized Jiittner-functions are used. 
However, in general, other momentum distributions could also be used. We neglect the 
melting of the prehadrons back into the deconfined phase. This equation assumes the 
escape of color-neutral prehadrons from the deconfined phase. 

In the MICOR model we determine the momentum spectra of the produced prehadrons 
generalizing a momentum dependent version of the rate equation found in eq. (|1|). We 
introduce an extended averaging on the r = th hadronization hypersurface: 

((-12)) = , (3) 

where 

/[...] = j dV,dV2d^p,£p2 [...] = (4) 

= J dV^dV2d^p^d'p2Q{E^)Q{E2){^ErE2)5{pl - ml)5{pl - ml) [...] (5) 

is the 12 dimensional phase space integration operator. Assuming energy-momentum 
conservation, the momentum of the outgoing particle is the sum of the incoming momenta: 
p = Pi + P2- Since we are interested in the momentum spectra of the outgoing particle, 
it is more useful to represent / in eq. (|^) as an integral of the relative four-momentum 
q = {pi — ^2)72 and the outgoing four-momentum p with s = p^: 

/[...] = Jd'pi,[...], (6) 
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Using the parametrization 

p^ = {m'rpchy, p'rp cos if, p'rpsmip, m^shy), (8) 
q^ = ±\/i^(9c, g^sinCcos(v9 + x), gsSinCsin((/) + x), g^cosC), (9) 
after some tedious algebra we obtain a simplified expression for the phase space integral: 



2vVl +1 , ... 

= / dV,dV2^\Y. I qsdqc I dcosC2E,E2^^^^^ [...] . (10) 
_^ J J ■'X-i-(cosC) 



^ e(g2) 



Here 



X-i-(cosC) = —acos^C + b±cosC — c±, (11) 
a = (E^-s)qI (12) 

PzQs, (13) 




Thus, from eq. (PH]), we derived the differential form of the averaging in eq. 

d{{(7Vl2)) ^ ip[flf2(TVl2] 

d'P /[/1/2] 



(15) 



The eq. (|T5D leads to a momentum dependent rate equation for the produced prehadrons 
and diquarks. 

In eq. (|l^) the prehadron is created off-shell, s = = m^^^^. When the prehadron 
leaves the deconfined region, it assumes to be an on-shell resonance by emitting or absorb- 
ing energy. For simplicity, we assume that its velocity distribution remains unchanged 
during this process. One can parametrize the transverse velocity as 



ch/i = mj./vs , (16) 
where = \Jpt + s. With this parametrization the differential momenta element is 

d'^p = sds chfidch^dydip, (17) 
and thus the four-velocity distribution is 

- fdss^^^^. (18) 



oh. jjL d ch. jj, dy dip J d^p 



If we substitute ch/x = rriT/mh into eq. ([T8|) and convert the velocity distribution 



into a four-momentum distribution, we obtain the on-shell momentum distribution for 



hadrons with mass nih and transverse mass rriT = vPt ~^ "^h- 



1 . (19) 



mTdrriTdy dip m\ J d^p 

This formula gives the momentum distribution of the primarily produced excited hadrons. 
From the expression in eq. (|l^) one can determine in one step the transverse momentum 
slope of the and in two steps the f2 (included the formation of an ss diquark). These 
results can be compared directly with the experimental data. For the description of all 
final state hadrons we would need to follow the time-evolution of the multicomponent 
hadronization and the decay of excited hadrons. This is beyond the scope of this paper 



and those calculations will be published in a following paper |jT8|. 



The deconfined matter is characterized by the hadronization temperature, Thadi and 
an initial transverse flow, f j,. We consider a large enough longitudinal extension for the 
deconfined matter: rjQ = ±1.8, where r]Q is the space-time rapidity. This parameter will 
not influence the transverse momentum spectra at y = 0. The constituent quark masses 
are chosen to be mg = 300 MeV and ms = 450 MeV. 

We vary the hadronization temperature in the region Th^d = 130 — 260 MeV and the 
initial transverse flow f ° = — 0.7. The obtained transverse spectra for the (p and the Q 
are fitted in the measured transverse momentum region 0.5 < rriT — Tn^ < 2.2 GeV for cf) 
and 0.3 < rriT — < 1.5 GeV for Q, following the procedure of the NA50 and WA97 
Collaboration, respectively. We compare the theoretical slopes to the experimental data. 
Fig. 1 displays the values of obtained from the measured and the calculated slopes of 
the (f) and the Q spectra. The area inside the solid contour line indicates < 3.67. 

One can see a very strong correlation between the hadronization temperature and the 
initial transverse flow, which can be characterized by the expression 

Thad + a ■ (4)2 = 0.25 GeV , (20) 

where a = 0.37 GeV. This correlation is indicated by the solid dark line in Fig. 1. 
Considering only the < 3.67 values, a large temperature region is allowed, namely 
Thad = 160 — 230 MeV (paired with the appropriate initial transverse flow f j.), for an 
initial condition for the hadronization. On the bottom part of Fig. 1 we show the details 
of our fit in an enhanced way, displaying the values of l/x^ < 3 on a lego plot. Here 
an 'excellent' agreement can be seen between the result of the MICOR model and the 
experimental data in the region Thad = 175 ± 15 MeV and = 0.46 ± 0.05. 

To avoid later confusion between the theoretical result of this paper and the experimen- 
tal data of NA50 and NA49 Collaboration |^, Fig. 2 displays our recent information 
on the transverse momentum spectra of the (p nieson. In the MICOR model, according 
to eq.(pOD we use Thad = 175 MeV and Vj' = 0.46 to calculate the theoretical transverse 
momentum spectra of the meson. 

The full squares show the NA49 data as reconstructed from Ref. The solid line 
indicates the T = 295 ±15 MeV slope with the parameterization 

dN/mTdrriT oc exp(— m^/T) (21) 

in the momentum region 0.02 < mx — < 1.5 GeV. 

The dashed line indicates the results of the NA50 Collaboration namely the 
T = 222 ± 6MeV slope with the parametrization 

dN/mTdrriT = C ■ ttit ■ Ki {niT/T) (22) 

in the momentum region 0.5 < mj- — < 2.2 GeV. The dotted line, which parallels the 
dashed one, indicates the same result as NA50, when using the expansion of the Bessel 
function 

dN/mTdrriT = C ■ ^ exp(— m^/T) , (23) 

where C = C ■ \JTtt/2. In this momentum region the difference between these two 
parametrization is negligible, and so the results of NA50 can be compared with the results 
of WA97 [Q, who used the parametrization of eq. (|23|) . 



Fig. 2 shows that the slopes of the NA49 Collaboration is valid for the transverse 
momentum region rriT — < 1 GeV and the slope of the NA50 Collaboration is valid for 
rriT — m^ > 0.3 GeV. The MICOR result was chosen to agree with the NA50 measurement, 
and it agrees with the experimental results of the NA49 in the region rriT — > 0.3 
GeV. In parallel, Fig. 2 demonstrates that the coalescence process creates a non-thermal 
momentum distribution, and the obtained momentum distribution for the (p can mimic 
a thermalized final state in large momentum region. (The fluctuations on the MICOR 
results are coming from the applied Monte-Carlo evaluation of the phase-space integral.) 
Fig. 2 shows that the MICOR model fails to reproduce the NA49 data in the momentum 
region rriT — fri^ < 0.3 GeV, which requires further improvements at small pr- However 
this disagreement will not influence our analysis and discussion at larger px- 

We apply the obtained correlation between the hadronization temperature and early 
transverse flow from eq. (|20| ) and recalculate the slopes of the transverse momentum spec- 
tra for Q and 0. Fig. 3. shows our results, displaying the "Data/Theory" ratios for the 
different Tgiope values (including the experimental error bars) as a function of hadroniza- 
tion temperature, Thad- This figure confirms the validity of the function used in eq.(^Dp. 

We also calculate the p/u spectra. In the MICOR model the p/u particle is formed in 
one step from the coalescence of a quark and an antiquark. The theoretical value of this 
slope parameter was fit in the measured transverse momentum region 1.5 < rriT < 3.2 
and compared to the experimental value T(p) = 219 ± 5 MeV. The bottom part of Fig. 3 
shows that the pairs of Thad and Vj^ parameters which satisfy eq. ( pOD by reproducing the 
slope of (p and Q, also reproduce the measured slope of p. This occurs especially in the 
region Thad = 175 ± 15 MeV. 

This result is very surprising. We could expect that the weakly interacting and Q will 
conserve their transverse momentum distribution from the initial hot phase (according to 
our assumption this is a strongly coupled deconfined phase). But the p particle is very 
strongly coupled to nucleons and pions [^], thus final state interactions should modify 
their spectra with an extra transverse boost. It was hardly expected, that the p meson 
could conserve any of his early transverse momentum distribution. 

With this surprising result, we collect the available data on the Tgiope extracted by the 
parametrization of eq. (^3]), similarly to p, and Q. We show these experimental data 
on Fig. 4. The data on Kg, A, A, S~, S and Q can be found in Ref. 0. The data on p 
and are from Ref. |T[]. For 7r+, tt^, K~^, , p^ and p~ we considered the NA44 data 



24] and using eq. (E3f) fit these data in the momentum region ttit — rrii > 0.3 GeV. We 



obtained T(7r+) = 150 MeV, T(7r-) = 145 MeV, T{K+) = 225 MeV, T{K-) = 216 MeV, 
T(p+) = 234 MeV, T(p") = 223 MeV. These data followed the tendency that the value 
of the Tsiope decreases with 25-30 MeV when using eq. (p3D as opposed to eq. (pTj), see e.g. 
Ref. 0. In Fig. 4, the dotted line was drawn across the p, and f2. 

Fig. 4 helps to estimate the efficiency of the final state interactions in the transverse 
momentum region ttit — rrti > 0.3 GeV. It reveals that the separation of the p"*", p~. A, 
A, and H from the weakly interacting and VL is much smaller than indicated in 
Refs. [|, 1^. Even more, according to our fit, the slope of p^ and p~ are very close to 
that of the weakly interacting meson. Thus the slope of the p meson and its restoration 
in the MICOR model is not surprising anymore. Fig. 4 indicates that in the transverse 
momentum region mT — Trii > 0.3 GeV the hadronic spectra may not suffer large modifica- 
tions because of secondary collisions. Pions could be an exemptions, since their separation 



from the other particles may indicate extra-long hadronic evolution, or their appearance 
from resonance decays. These results demand further investigation into the properties of 
the hadronic spectra at rriT — rrii < 0.3 GeV. 

In this paper, we analyzed the recent experimental results on the production of the 
meson and the Q baryon, using the MICOR model. We have found a strong correlation 
between the hadronization temperature and the transverse flow of the early deconfined 
state, where the and Q were created from quark- antiquark plasma. We have found that 
Thad = 175 ±15 MeV and v° = 0.46 ±0.05 is favoured by data. In this temperature region 
the MICOR model also reproduced the measured p/u slope. This result indicates that 
final state interactions modify very weakly the transverse momentum spectra of the p/u 
meson in the transverse momentum region 1.5 < ttit < 3.2 GeV. In the large mj- region 
the final state interactions may not modify the transverse slopes created in an early state. 
Further calculations are needed to confirm if all hadronic slopes can be reproduced by the 
MICOR model which assumes quark coalescence. A positive result may be interpreted as 
the formation of a deconfined phase (namely a massive quark-antiquark plasma) in the 
early stage of the PbPb collision at CERN SPS. 
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Fig. 1. The x^"fit of the calculated T^ff slopes for the (f) meson and the Q baryon. 
The area inside the solid contour line indicates < 3.67, The solid dark line shows 
the fit Thad + a ■ {vt)^ = 0-25 GeV, where a = 0.37 GeV. On the bottom figure we 
show our fit in an enhanced way, displaying the values of 1/%^ < 3 on a lego plot. 
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Fig. 2. The transverse momentum spectra of the meson measured by the NA49 
Collaboration in Pb+Pb collision at 158 AGeV (full squares) 0. The solid dark 
line indicates the effective slope T = 295 ± 15 MeV published by NA49. The dashed 
line (and the overlapping dotted line) indicates the measured slope of the NA50 
Collaboration |l[] in the transverse momentum region 0.5 < ttit — < 2.2 (for 
explanation see text). The thin full curve with fluctuations is the MICOR result for 
production calculated at Thad = 175 MeV and Vj^ = 0.46. 
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Fig. 3. The "Data/ Theory" values for recalculated Tsiope{fi) and Tsiope{<P) as a func- 
tion of hadronization temperature T^ad with transverse flow satisfying eg. (PO]) . 
We displayed the same ratio for the p particle calculated in the MICOR model 
(bottom part). Vertical dotted line indicates T^ad = 175 MeV. 
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Fig. 4. Experimental liadronic slopes of the transverse momentum spectra in the 
PbPb collision at 158 AGeV energy from WA97 (squares), NA50 (dots) and 
NA44 Collaboration (triangulars) . The results of WA97 and NA50 were fitted 
originally by eq.(^) and we fit the NA44 data on tt^, K"^, p"*" and p~ ||2^ in the 
same way in the momentum region rriT — frii > 0.3 GeV. 



